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A solar-powered high temperature differential Stirling engine has been considered for optimization with 
multiple criteria. A mathematical model based on the finite-time thermodynamics has been developed so 
that the output power and thermal efficiency and the rate of entropy generation of the solar Stirling sys¬ 
tem with finite rate of heat transfer, regenerative heat loss, conductive thermal bridging loss and finite 
regeneration process time are obtained. Furthermore, imperfect performance of the dish collector and 
convective/radiative heat transfer mechanisms at the hot end as well as the convective heat transfer at 
the heat sink of the engine are considered in the developed model. Three objective functions including 
the output power and overall thermal efficiency have been considered simultaneously for maximization 
and the rate of entropy generation of the Stirling engine are minimized at the same time. Multi-objective 
evolutionary algorithms (MOEAs) based on NSGA-II algorithm has been employed while the Effective¬ 
ness’s of regenerator, the Effectiveness’s of the low temperature heat exchanger, the Effectiveness’s of 
the high temperature heat exchanger, heat capacitance rate of the heat sink, heat capacitance rate of 
the heat source, temperatures of the working fluid in the high temperature isothermal process and tem¬ 
peratures of the working fluid in the low temperature isothermal process are considered as decision vari¬ 
ables. Pareto optimal frontier has been obtained and a final optimal solution has been selected using 
various decision-making approaches including the fuzzy Bellman-Zadeh, LINMAP and TOPSIS methods. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine is a simple type of external-combustion en¬ 
gine that uses a compressible fluid as a working fluid. The Stirling 
engine can theoretically be a very efficient engine to convert heat 
into mechanical work at Carnot efficiency. The thermal limit for 
the operation of a Stirling engine depends on the material used 
for its construction. In most instances, the engines operate with a 
heater and cooler temperature of 923 and 338 K, respectively [1]. 
Engine efficiency ranges from about 30% to 40% resulting from a 
typical temperature range of 923-1073 K, and normal operating 
speed range from 2000 to 4000 rpm [2]. 

Urieli and Berchowitz [3], Reader and Hooper [4] and Har¬ 
greaves [5] among others all provide isothermal models similar 
to Schmidt’s original work. Carlson et al. [6] developed an ideal 
model with non-isothermal heat exchange. Models of this type 
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suggest an improvement on the isothermal model as they elimi¬ 
nate the necessity for infinite heat transfer and impractically slow 
engine speed associated with isothermal working spaces. Urieli 
and Kushnir [7] showed that this analysis can be utilized in order 
to evaluate the various practical effects of non-ideal regenerators, 
heat exchangers, including heat transfer and pressure losses. Mar- 
taj et al. [8] presented a thermodynamic analysis of a low temper¬ 
ature Stirling engine at steady state operation, and energy, entropy 
and exergy balances were presented at each main element of the 
engine. The major goals of the Stirling engine designers can be 
mentioned in three categories: maximum efficiency; maximum 
power; minimum costs. Markman et al. [9] studied the thermal- 
flux and mechanical-power losses of a 200 W beta-configuration 
of the Stirling engine to optimize and increase the engine 
efficiency. 

A miniature Stirling engine with 4 W output power in 0.1 MPa 
pressure and 900 rpm rotation speed was studied and made by 
Kagawa et al. [10]. Brandhorst and Chapman [11] developed a 
5 kW engine for use as a power generator in space applications. 
Ataer [12] employed the Lagrangian method to the analysis of 
regenerators of Stirling cycle machines. The governing equations 















M.H. Ahmadi et al.) Energy Conversion and Management 75 (2013) 438-445 


439 


Nomenclature 



A 

area (m 2 ) 

1,2,3,4 

state points 

C 

heat capacitance rate (W K _1 ) 

2 

outlet 

Cv 

specific heat capacity (J mol -1 K _1 ) 

app 

collector aperture 

h 

heat transfer coefficient (W m -2 K” 1 ) 

ave 

average 

I 

direct solar flux intensity (W m -2 ) 

c 

cold side 

K 0 

heat leak coefficient (W K 1 ) 

dis 

distance 

M 

proportionality constant (-) 

H 

heat source 

n 

number of mole 

h 

hot side 

N 

number of heat transfer units (-) 

L 

cold side/heat sink 

P 

power output (W) 

R 

regenerator 

Q 

heat (J) 

rec 

absorber 

R 

the gas constant (J mol -1 K _1 ) 



S 

entropy (J/K) 

Greek 


T 

temperature (I<) 

p 

thermal efficiency [-] 

t 

time (s) 

8 

effectiveness and emissivity factor [-] 

U 

overall heat transfer coefficient (W K _1 m -2 ) 

A 

ratio of volume during the regenerative processes [-] 

V 

volume of the working fluid (m 3 ) 

3 

Stefan’s constant [W m -2 K -4 ] 

w 

output work (j) 

G 

entropy production [WK' 1 ] 

Subscripts 



0 

ambient 



1 

inlet 




of the regenerator are derived in terms of the displacement of the 
displacer, so that time does not appear in the equations. The equa¬ 
tions, which include pressure fluctuations owing to flow reversals 
and longitudinal conduction, are solved numerically by a digital 
computer using a finite difference method. Nakajima et al. [13] 
developed a 10 g micro-Stirling engine with an approximately 
0.05 cm 3 piston swept volume. An engine output power of 
10 mW at 10 Hz was reported. The problems of scaling down were 
discussed. 

Aramtummaphon [14] examined an open cycle Stirling engines 
by using steam heated from producer gas. The first engine gener¬ 
ated an indicated power of about 1.36 kW at a maximum speed 
of 950 rpm, while the second engine, improved from the first 
one, produced an indicated power of about 2.92 kW at a maximum 
speed of 2200 rpm. Fukui et al. [15] has designed and constructed a 
micro-engine, and its experimental examination was performed; 
however, the micro-engine performance cannot be scaled to prac¬ 
tical engine. Experimental determination of the effect of operating 
variables on Stirling engine is complicated and time consuming. 
Iwamoto et al. [16] comprised low temperature and high temper¬ 
ature Stirling engine efficiency with others. They determination 
exhibited that efficiency of LTD Stirling engine is around 50% of 
Carnot efficiency with same situation. Wu et al. [17] showed the 
effects of heat transfer, regeneration time, and imperfect regener¬ 
ation on the performance of the irreversible Stirling engine cycle. 
Erbay and Yavuz [18] analyzed the real Stirling heat engine for 
maximum power output conditions using polytropic processes. 
They also determined the efficiency and compression ratio at max¬ 
imum power density and ascertained the thermal design bounds. 
Ahmadi and Hosseinzade [19] investigated of Solar Collector 
Design Parameters Effect onto Solar Stirling Engine Efficiency. 
Ahmadi et al. [20] developed intelligent approach to figure power 
of solar Stirling heat engine by implementation of evolutionary 
algorithm. 

Yaqi et al. and Sharma et al. developed a mathematical model 
for the overall thermal efficiency of solar powered high tempera¬ 
ture differential dish Stirling engine with finite heat transfer and 
irreversibility of regenerator and optimized the absorber 
temperature and corresponding thermal efficiency [21,22]. Tlili 


investigated effects of regenerating effectiveness and heat capaci¬ 
tance rate of external fluids in heat source/sink at maximum power 
and efficiency [23]. Kaushik et al. studied the effects of irreversibil¬ 
ities of regeneration and heat transfer of heat/sink sources [24-27]. 
Solution of the multi-objective optimization problems is an extre¬ 
mely difficult goal which requires the simultaneous satisfaction of 
a number of different and even conflicting objectives. Evolutionary 
algorithms (EA) were initially extended and employed during the 
mid-eighties in an attempt to stochastically solve problems of this 
generic class [28]. A reasonable solution to a multi-objective prob¬ 
lem is to investigate a set of solutions, each of which satisfies the 
objectives at an acceptable level without being dominated by any 
other solution [29]. Multi-objective optimization problems in gen¬ 
eral show a possibly uncountable set of solutions namely as Pareto 
frontier, whose evaluated vectors represent the best possible 
trade-offs in the objective function space. In this term, multi-objec¬ 
tive optimization of different themodynamic and energy systems 
have been paid attention by researchers nowadays [30-37]. In this 
work, by implementing multi-objective optimization algorithms 
output power and Stirling engine’s thermal efficiency are maxi¬ 
mized and rate of entropy generation of the Solar-Dish Stirling en¬ 
gine minimized. Also final solutions of different multi-objective 
optimization were compared with Ref. [24] data. At end it should 
be mentioned that, error analysis were implemented to figure 
robustness and precision of final solutions of different decision 
making approaches out. 


2. System description 

In Solar-Dish Stirling systems, mirrors of the parabolic shaped 
concentrator focuses the sun light on the focal point of the concen¬ 
trator where the hot end of the Stirling engine is located. Therefore, 
the solar energy with a relatively high temperature is transferred 
to the hot side heat exchanger of the Stirling engine. Fig. 1 illus¬ 
trates a schematic for a Solar-Dish Stirling engine connected to a 
solar dish concentrator. The Solar-Dish is equipped with a sun 
tracker which tracks the sun in order to have maximum solar en¬ 
ergy transfer to the engine when sun moves during days. Hence, 
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Fig. 1. Schematic of a solar Stirling engine [32]. 
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S 




Fig. 3. T-S diagram of a Stirling engine cycle. 


solar energy is absorbed and transferred to the working fluid in the 
hot space of engine. 

Fig. 2 is a schematic diagram of a Stirling heat engine cycle with 
finite-time heat transfer and regenerative heat losses. As shown in 
Fig. 3, Stirling cycle consists of four processes. Process 1-2 is an iso¬ 
thermal process, in which the working fluid after compressing at 
constant temperature, T c and rejected heat to the heat sink at 
low temperature T Ll . Therefore, the temperature of heat sink is in¬ 
creased to Tl 2 . Then the working fluid crosses over the regenerator 
and warms up to T h in an isochoric process 2-3. In process 3-4, the 
working fluid is expanded in a constant temperature T h process and 
receives heat from the heat source in which its temperature is re¬ 
duced from Th, into W Last process (4-1), is an isochoric cooling 
process, where the regenerator absorbs heat from the working 
fluid. 

3. Thermodynamic analysis of the system 

Actual useful heat gain of the dish collector, considering con¬ 
duction, convection and radiation losses is given by [20-22] as 
follows: 

Qu = Mo - Aeclh(T H „ e - To) + eS(T 4 m - T <)] (1) 

where I is the direct solar flux intensity, A app is the collector aper¬ 
ture area, rj 0 is the collector optical efficiency, A rec is the absorber 
area, h is conduction/convection heat transfer coefficient, T Have is 
the average absorber temperature, T 0 is the ambient temperature, 
s is emissivity factor of the collector, 3 is the Stefan’s constant. 

Thermal efficiency r] s of the dish collector is obtained as 
[ 20 - 22 ]: 

% = ir~ = n 0 - ^ - To) + eS(T 4 Hm - Tg)] (2) 

jCi app * C 


3 A. Regenerative heat losses in the regenerator 


It is important to mention that there also exists a finite heat 
transfer in the regenerative heat transfer (Qr) which is given by 
[ 20 - 22 ] 

Q r = nC„s R (T h - T c ) (3) 


Where A Qr is the heat loss during the two regenerative processes in 
the cycle. The following relationship is obtained [20-22]: 

AQ r = nC,(l -s R )(T h -T c ) (4) 


n is the mass of the working fluid in mole, C v is the specific heat 
capacity of the working fluid in the regenerative processes in terms 
of mole, £ r is the effectiveness of regenerator, T h and T c are the work¬ 
ing fluid temperatures in the hot space and cold space, respectively. 

Owing to the influence of irreversibility of the finite-rate heat 
transfer, the time of the regenerative processes is not negligible 
in comparison to that of the two isothermal processes [20]. In or¬ 
der to calculate the time of the regenerative processes, one as¬ 
sumes that the temperature of the working fluid in the 
regenerative processes as a function of time is given by [20]: 


dT 

dt 


±Mi 



where M is the proportionality constant which is independent of 
the temperature difference and dependent only on the property of 
the regenerative material, called regenerative time constant and 
the ± sign belong to the heating (i = 1) and cooling (i = 2) processes 
respectively [20]. 
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3.2. The amounts of heat released by the heat source and absorbed by 
the heat sink 


Considering the cyclic period of the Stirling engine, the output 
power, the thermal efficiency and entropy production of the engine 
are given by: 


W _ Q H ~ Ql 

t t 



The heat released between heat source and working fluid (Qh), 
the heat absorbed between the working fluid and the heat sink 
(Qc) is obtained as follows 


Q h = nRT„LrU + nC v ( 1 - s R )(T h - T c ) 

(8) 

Qc = nRTcLn/l + nC v ( 1 - e R )(T h - T c ) 

In the other hand: 

(9) 

Qh = Ch£h(^H! - + £ C h £h(Thi ~ T^j th 

(10) 

Qc = Cl£l(T c - T’lJG 

(11) 


dt = 


Qh-Ql 


1 

G = - 


Qh 
Ql q 


H 


t \Ti T 


-ave 


H, 


( 20 ) 


( 21 ) 


ave, 


Substituting Eqs. (3)-(14) into Eqs. (19) and (20) we have, 

nR(T h -T c )Ln2 


P = 


nRT h LnA+nC v (t-£ R )(T h -T c 


Chs h (t h1 -r h )uc H s H (r^-rA 


nRT c lnX+nC v A-£ R )(T h -T c ) , ( 1 , 1 \ (t T\ 
C L e L (T c -T h ) +\M,+M 2 )Vh 1 c) 


( 22 ) 


where C H and C L are the heat capacitance rate of external fluids in 
the heat source and heat sink, respectively. 


s H = 1 - e Nh 

(12) 

8 L = 1 - e~ Nl 

(13) 


where s H and s L are the Effectiveness’s of the high and low 
temperature heat exchangers, respectively. And N L = U L A L /C L , 
N h = U h A h /C h are the cyclic period. Using Eqs. (3)-( 11), we get that 
the cyclic period t is: 

t _ nRT h Ln2 + nC„(l -e R )(T h - T c ) 

~ C H s H (T Hl - T h ) + ZC H e H (Tl - T 4 ,) 

nRT c LnA + nC„(l - E R )(T h - T c ) (\ IV 
C l e l (T c -T L i ) + {MiM 2 J [h 

- Tc) (14) 

3.3. The conductive thermal bridging losses from heat source to the 
heat sink 


This value is proportional to the average temperature difference 
of the heat source and heat sink and the cycle time, it is obtained as 
follows: 


Qo — Kq {TH ave Tl ave ) t cycle 

(15) 

T _ Thi + Th 2 

■* Have 2 

(15a) 

T _ Tl x + Tl 2 

tave 2 

(15b) 

We have [24-27] 


Th 2 = (1 - £h)Th x + SnTh 

(15c) 

Tl 2 = (1 — £l)Tl x + £lT c 

Thus, using Eqs. (15a)-(15d), we have 

(15d) 

Qo = 2 [(^ £ h)Th x -( 2 £l)Tl x + (GhT h £lT c )] t cyc \ e 

(16) 


The net heat released from the heat source (Qh) and the net heat 
absorbed by the heat sink (Qi) are obtained as follows: 


Qh = Qh + Qo (17) 

(18) 


_ nR(T h -T c ) Ln2 

^ nRT/,Lnl + nCV(l — £r)(Ti 2 — T c )+^ l [(2 — Sh)T Hi —(2 — Sl)T Li + (s^T \ c)]t C ycie 

(23) 

The maximum thermal efficiency of the entire Solar-Dish Stir¬ 
ling engine is product of the thermal efficiency of the collector 
and the optimal thermal efficiency of the Stirling engine [22]. 
Namely: 

Vm = fls'lt ( 24 ) 

4. Multi-objective optimization with evolutionary algorithms 

4.1. Optimization via EA 

In this study, the Pareto frontier is obtained using the Genetic 
Algorithm (GA) known as a branch of evolutionary algorithm. Ge¬ 
netic algorithms were developed by John Holland in the 1960s as a 
means of importing the mechanisms of natural adaptation into 
computer algorithms and numerical optimization [31]. They are 
implemented as a computer simulation in which a population of 
abstract representations (called chromosomes or the genotype of 
the genome) of candidate solutions (called individuals, creatures, 
or pheno-types) to an optimization problem evolves toward better 
solutions. The evolution normally starts from a population of ran¬ 
domly generated individuals and happens in generations. In each 
generation, the fitness of every individual in the population is eval¬ 
uated; multiple individuals are stochastically selected from the 
current population (based on their fitness), and modified (recom¬ 
bined and possibly randomly mutated) to form a new population. 
The new population is then used in the next iteration of the algo¬ 
rithm. Commonly, the algorithm terminates when either a maxi¬ 
mum number of generations have been produced, or a 
satisfactory fitness level has been reached for the population. If 
the algorithm has terminated due to a maximum number of gener¬ 
ations, a satisfactory solution may or may not have been reached. 
In genetic algorithms, a candidate solution to a problem is typically 
called a chromosome, and the evolutionary viability of each chro¬ 
mosome is given by a fitness function. This method is a powerful 
optimization tool for nonlinear problems [29,32]. 

Also, Multi-objective evolutionary algorithms (MOEAs) have 
been developed over the past decade by numerous tests on com¬ 
plex mathematical problems and on real-world engineering prob¬ 
lems and have shown that they can eliminate the difficulties of 
classical methods [29,32]. The structure of the MOEA used in the 
present study is illustrated in Fig. 4 [31 ]. The real values of decision 
variables are used instead of their binary coded. 


Ql — Q-c + Q-o 
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Fig. 4. Scheme for the multi-objective evolutionary algorithm used in the present 
study [31-33]. 

4.2. Non-dominated sorting genetic algorithm, NSGA-II 

NSGA-II algorithm is utilized here as a multi-objective optimi¬ 
zation method in order to find the Pareto frontier using the genetic 
algorithm. In this algorithm solutions are categorized based on the 
Pareto concept and sorting non-dominated solutions into non- 
dominated layers as schematically shown in Fig. 5. In other words, 
if N p is the number of population, it is categorized into N L layers in 
which intersection of each two arbitrary selected layers is empty 
set and union of all layers is N p set. 

Elitism or virtual fitness of each individual is equal to its layer. 
Parent selection for cross over operation is performed based on the 
tournament selection between two random selected layers. There¬ 
fore, individual located on the layer 1, have more chance to be se¬ 
lected for the next generation. Uniform distribution of solutions 
along layers is controlled with introducing of crowding distance in¬ 
dex for each solution. This index is defined as a ratio of subtraction 
of objective functions for two neighbor solutions around the cur¬ 
rent solution to the subtraction of the maximum and minimum 
values of that objective. Therefore, for /<th objective of jth solution, 
we have, 

j fkj-l ~fkj +1 nn 

IdiSjJi = ~r - —r - (/~>J 

Jk, max Jk, min 

For boundary solutions (solutions with smallest and largest 
function values) are assigned an infinite distance index. The overall 
crowding distance value is calculated as the sum of individual dis¬ 
tance values corresponding to each objective. 

M 

Idisj — 'y^hlis.j.k (26) 

k= 1 



Fig. 5. Schematic of solution layering in NSGA-II algorithm [33]. 



Fig. 6. Schematic of distance indexing of individuals in NSGA-II algorithm [33]. 

where M is the number of objectives and j is the individual index. 
Fig. 6 shows schematic of evaluation of distance index. In this algo¬ 
rithm, two parameters are calculated for each solution: 

(1) Dominant number, NL, that is the number of solutions that 
dominate the current solution. Detail of dominancy concept 
and definition is explained well in multi-objective optimiza¬ 
tion in Refs. [28,30]. Dominant number, for non-dominated 
individuals of the current population is zero, therefore, these 
solutions are placed in layer 1. Non-dominated solutions for 
a set of the individuals excluding the layer 1 members are 
placed in layer 2. For an M objectives problem with N popu¬ 
lations, the number of comparisons is MN2. This procedure 
continued so that all individuals are accommodated in their 
appropriate layers. In the continuing, irank index for each 
individual is assigned as its dominant number (layer num¬ 
ber), NL. 

(2) Crowded comparison operator, n^, defined as follow: 

A Biff ran k A < rank B )Or : (frank,, = rank B )and \ disA > I dM ) (27) 

It means that for two individuals with differing non-dominated 
ranks (different layers), the solution with the lower rank (layer) is 
preferred. Otherwise, for two solutions of the same layer, the solu¬ 
tion placed in the region with lower concentration of solutions is se¬ 
lected. This is the “concept of dominated sorting algorithm.” 

4.3. Objective functions, decision variables and constraints 

Three key objective functions for optimization are the system’s 
entropy production (should be minimized), the output power 
(should be maximized) and the thermal efficiency of Dish-Stirling 
engine system’s (should be maximized), denoted by Eqs. (21), 
(22) and (24), respectively. 

In this paper, seven decision variables have been considered as 
follow: 

s R is the Effectiveness’s of regenerator; s L the Effectiveness’s of 
the low temperature heat exchanger; s H the Effectiveness’s of the 
high temperature heat exchanger; C L the heat capacitance rate of 
the heat sink; C H the heat capacitance rate of the heat source; T h 
the temperatures of the working fluid in the high temperature iso¬ 
thermal process 3-4; and T c the temperatures of the working fluid 
in the low temperature isothermal process 1-2. 

The objective functions with respect to following constraints 


have been solved: 


0.4 ^ s R ^ 0.9 

(28) 

0.4 ^ s H ^ 0.8 

(29) 
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Table 1 

Specified genetic algorithm options for optimization problem in this study. 


Specified options for GA 

Value 

Population type 

Duble vector 

Population size 

400 

Selection process 

Tournament 

Tournament size 

2 

Mutation 

Constraint dependent 

Maximum number of generations 

1000 


120 


100 




80 

60 


40 


20 

0.37 



Fig. 7. Pareto optimal frontier in objectives’ space. 


0.4 < s L < 0.8 (30) 

300 1800 (31) 

300 < C L < 1800 (32) 

800 < T h < 1000 K (33) 

400 ^ T c ^ 510 K (34) 


To find the optimum design parameters of the system, based on ge¬ 
netic algorithm technique a simulation program is developed in 
Matlab software. Genetic algorithm options for optimization prob¬ 
lem are listed in Table 1. 

4.4. Decision-making in the multi-objective optimization 

In multi-objective optimization a process of decision-making for 
selection of the final optimal solution from available solutions is re¬ 
quired. A number of several methods for decision-making process 
in decision problems exist. And these methods can be employed 
in decision-making for selection of a final optimal solution from 
the Pareto frontier which is obtained by optimization. In this paper 
most recognized and common type of decision-making processes 
including the fuzzy Bellman-Zadeh, LINMAP and TOPSIS method 
is used in parallel and the final optimal solution has been decided 



Fig. 8. Pareto optimal frontier in objectives’ space (rj m - P). 

based on engineering experience and criteria among solutions 
which suggested by these three methods. For more information 
about these methods see Refs. [38,39]. 


5. Result and discussion 

The output power and thermal efficiency of Solar Dish-Stirling 
system are maximized simultaneously while the rate of entropy 
generation of the Stirling engine are minimized simultaneously 
using the multi-objective optimizing method which works based 
on the NSGA-II algorithm. In this regard, optimization is performed 
with objective functions that are expressed by Eqs. (21), (22) and 
(24), also constraints are expressed with Eqs. (28)-(34). Design 
parameters (decision variables) of optimization are the Effective¬ 
ness’s of regenerator, Effectiveness’s of the low temperature heat 
exchanger, Effectiveness’s of the high temperature heat exchanger, 
heat capacitance rate of the heat sink, heat capacitance rate of the 
heat source, temperatures of the working fluid in the high temper¬ 
ature isothermal process 3-4, temperatures of the working fluid in 
the low temperature isothermal process 1-2. In order to have con¬ 
sistency with previous works, specifications of the Stirling engine 
are considered as follows [22,24], 

/ = 1000 W/m 2 , C = 1300, e = 0.9, t] 0 = 0.9, K 0 = 2.5 W/K, 
n = 1, C v = 15 J mol -1 K" 1 , R = 4.3Jmor 1 K"\ r H ,= 1300I<, 

T t , = 290 K, T 0 = 288 I(,{ = 2x 10~ 10 , h = 20 W rrr 2 K ’, 

2 = 2, 1 /Mi + 1 /M 2 = 2 x 10’ 5 s/K, S = 5.67 x 10“ 8 W itt 2 IC 4 . 

Fig. 7 illustrates the Pareto frontier in the proposed objectives’ 
space obtained using the multi-objective evolutionary algorithm 
(NSGA-II). Three final solutions selected by the Fuzzy Bellman-Za¬ 
deh, LINMAP and TOPSIS decision-makers are indicated separately 
in Fig. 7. 

Table 2 compares optimal results obtained in this paper using 
three decision making methods with corresponding results ob¬ 
tained in Ref. [24]. Where Table 1 demonstrated obtained results 


Table 2 

Decision making of multi-objective optimal solutions. 


Optimal solution 

Decision making method 

Decision variables 





Objective functions 


£r 



c„ 

G 

T h (I<) 

T c (K) 

dm (%) 

P (1<W) 

a (WIC 1 ) 

This work 

Bellman-Zadeh 

0.9 

0.8 

0.8 

1410 

774 

996.7 

400.5 

36.33 

37.37 

70 


LINMAP 

0.9 

0.8 

0.8 

1413 

823 

996.7 

400.5 

36.36 

38.86 

69.5 


TOPSIS 

0.9 

0.8 

0.8 

1424 

1252 

996.9 

400.5 

36.56 

49.64 

87.47 

Ref. [24] 

- 

0.9 

0.8 

0.8 

1000 

1000 

962.2 

462.2 

34.2 

69.97 

- 
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Fig. 9. Pareto optimal frontier in objectives’ space (er - rj m ). 



Fig. 10. Pareto optimal frontier in objectives’ space (cr - P). 


from different decision making approaches for Solar Dish-Stirling 
engine and results exhibited by Kaushik [24]. 

In Fig. 8 plot of the Pareto optimal frontier for objective func¬ 
tions of output power and Solar Dish-Stirling system’s thermal effi¬ 
ciency which final solutions determined by different decision 
making approaches are highlighted. As can be seen in Fig. 8, the 
output power decreases with increasing the thermal efficiency of 
the engine. To assist the optimal design of the Solar Dish-Stirling 
system, fitted curve (rj m = -3.242 x 10 -10 P 2 -3.458 x 10~ 5 

P+ 0.5569), which is valid in the range of 59 kW< P< 60.5 kW, is 
derived for the Pareto optimal points curve. 

A previous analogue was carried out to illustrate Pareto 
optimal frontier for objective functions of Solar Dish-Stirling 
system’s thermal efficiency and the rate of entropy generation 
of the Stirling engine. Output results from this analogue are 
demonstrated in Fig. 9 in which by increasing the thermal effi¬ 
ciency, rate of entropy generation increases. Also, a curve is fitted 


(cr = 4.377 x 10 14 85 + 19.87) to all the points obtained by 

multi-objective optimization. In Fig. 10 Pareto optimal frontier 
for objective functions of output power and the rate of entropy 
generation of the Stirling engine is depicted. Form Fig. 10 can 
be observed that, the rate of entropy generation of the engine in¬ 
creases when the output power increases. Moreover, a curve is 
fitted (<r = 6.051 x 1(T 5 P 2 + 0.001305 P+13.23) to the points ob¬ 
tained by multi-objective optimization. 

Mean absolute percentage error was carried out in error analy¬ 
sis of implemented methods. To assess this goal, 30 runs of each 
approach were performed to get final solution by Bellman-Zadeh, 
LINMAP and TOSIS decision making methods. Then, magnitude of 
each objective function such as output power, Solar Dish-Stirling 
system’s thermal efficiency and rate of entropy generation of the 
Stirling engine were compared with correspond experimental val¬ 
ues. First row of Table 3 demonstrates maximum absolute percent¬ 
age error (MAAE) of three decision making methods. Moreover, 
second row of Table 3 reports mean absolute percentage error 
(MAPE) of mentioned methods. 

6. Conclusion 

In this study, finite-time thermodynamics has been applied to 
determine output power, thermal efficiency and rate of entropy 
generation of the Solar Dish-Stirling engine. Factors such as dish 
collector performance, conductive and radiative heat transfer 
mechanism at heat source and heat sink of the engine are involved 
in the analysis. The output power, thermal efficiency and rate of 
entropy generation of the engine were considered simultaneously 
for multi-objective optimization the Effectiveness’s of regenerator 
(£k), the Effectiveness’s of the low temperature heat exchanger 
(e L ), the Effectiveness’s of the high temperature heat exchanger 
(e H ), heat capacitance rate of the heat sink (C L ), heat capacitance 
rate of the heat source (C H ), temperatures of the working fluid in 
the high temperature isothermal process (T H ) and temperatures 
of the working fluid in the low temperature isothermal process 
(T l ) were considered as design parameters. Multi objective evolu¬ 
tionary algorithm is considered based on NSGA-II algorithm and 
the Pareto optimal frontier in objectives space was obtained. A fi¬ 
nal optimal solution was selected from solutions of the Pareto fron¬ 
tier using three decision making methods including LINMAP, 
TOPSIS and fuzzy methods. 

It was shown that multi objective-multi variable approach leads 
to more desirable design of engine if it is compared with traditional 
single objective-single variable approach. The priority of the multi 
objective-multi variable approach over the traditional single objec¬ 
tive-single variable was proven as in all cases significant lower 
deviation indexes was obtained for multi objective-multi variable 
approach. Further in this case it was observed that although the 
TOPSIS decision-making method yields a final solution with more 
desirable output power and thermal efficiency, the rate of entropy 
generation is increased. Error analysis was performed based on the 
MAPE method and it was shown that the average error of solutions 
obtained through three decision making methods were 1.3%, 4.4% 
and 3.5% for the output power, thermal efficiency and rate of entro¬ 
py generation, respectively. The maximum error of solutions ob¬ 
tained by three decision making methods were 2.5%, 8.4% and 


Table 3 

Error analysis based on the mean absolute percent error (MAPE) method. 


Decision making method 

Bellman-Zadeh 


LINMAP 



TOPSIS 



Objectives 

i] m (%) P (kW) 

d ( wr 1 ) 

dm (%) 

P (kW) 

(7 (W) 

dm (%) 

P (kW) 

(7 (WK" 1 ) 

Max Error (%) 

2.5 6.8 

6.2 

2.5 

6.5 

6.8 

2.4 

8.4 

6.5 

Average Error (%) 

1.3 2.4 

3.4 

1.2 

2.7 

3.0 

1.3 

4.4 

3.5 
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6.8% for the output power, thermal efficiency and rate of entropy 
generation, respectively. 
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